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FOREWORD 

This  report  covers  theoretical  work  on  supersonic  radial  outflow 
systems  for  investigating  radial  flow  diffusers,  carried  out  under 
Work  Unit  No.  2307N435  of  Project  No.  2307  in  the  Thermomechanics  Branch, 
Aeromechanics  Division,  Flight  Dynamics  Laboratory,  from  November  1977 
to  November  1978. 


iii 


AFWAL- TR- 80- 3028 


TABLE  OF  CONTENTS 

SECTION  PAGE 

I  INTRODUCTION  1 

II  RING  NOZZLE  DESIGN  2 

1.  Description  of  the  Nozzle  System  2 

2.  Derivation  of  the  Ring  Nozzle  Design  Equations  3 

3.  Numerical  Evaluation  of  the  Design  Equations  7 

4.  Discussion  of  the  Ring  Nozzle  Design  Chart  7 

5.  Design  Examples  9 

III  Ring  Nozzle  Performance  Analysis  14 

1.  General  Considerations  14 

2.  Derivation  of  the  Performance  Equations  15 

3.  The  Wall  Force  Parameter  20 

4.  Flow  Divergence  22 

IV  WALL  FRICTION  LOSSES  IN  AN  EXPANSION  NOZZLE  24 

1.  General  Considerations  24 

2.  Characteristic  Flow  Conditions  for  Wall  Friction 

Losses  24 

3.  Non-Circular  Expansion  Nozzle  26 

4.  Wall  Friction  Coefficient  for  an  Expansion  Nozzle  28 

V  Theoretical  Performance  of  a  Radial  Outflow  System  32 

1.  General  Aspects  32 

2.  Performance  Calculations  32 

3.  Accuracy  of  the  Performance  Calculations  34 

v 


f 


AFWAL-TR-80- 3028 


TABLE  OF  CONTENTS  (CONCLUDED) 

SECTION  PAGE 

VI  RADIAL  DIFFUSER  DESIGN  CONSIDERATIONS  36 

1.  The  Multi-Channel  Radial  Diffuser  36 

2.  Specific  Diffuser  Design  Considerations  37 

3.  The  Straight  Flow  Supersonic  Diffuser  38 

VII  CONCLUSIONS  44 


vi 


AFWAL- TR- 80- 3028 


LIST  OF  ILLUSTRATIONS 

FIGURE  PAGE 

1  Scheme  of  a  Ring  Nozzle  Expansion  System  50 

2  Ring  Nozzle  Design  Chart  for  Helium  51 

3  Ring  Nozzle  Design  Chart  for  Air  52 

4  Ring  Nozzle  Design  Chart  for  CO^  53 

5  Ring  Nozzle  Design  Chart  for  Freon  12  54 

6  Design  of  a  Ring  Nozzle  Expansion  System  for  M  =  2.3  55 

7  Scheme  of  a  Radial  Outflow  System  56 

8  Deflection  of  the  Radial  Flow  into  Individual 

Diffuser  Channel  57 

9  Distribution  of  the  Wall  Friction  Forces  Along  a 
Supersonic  Expansion  Nozzle  for  a  Constant  Wall 

Friction  Coefficient  58 

10  Diffuser  Plate  Design  for  a  Supersonic  Multi-Channel 

Radial  Diffuser  59 

11  Wall  Contour  and  Cross-Sectional  Area  for  Two  Different 

Radial  Diffuser  Inlet  Geometries  60 

12  Inlet  and  Exit  E-Function  According  to  Equation  80  for  a 

Diverging  Diffuser  Channel  61 

13  Inlet  and  Exit  E-Function  According  to  Equation  80  for  a 

Constant  Area  Diffuser  Channel  62 

14  Inlet  and  Exit  E-Function  According  to  Equation  80  for  a 

Converging  Diffuser  Channel  63 


vii 


/ 


AFWAL-TR-80- 3028 


SECTION  I 
INTRODUCTION 

A  preferred  flow  arrangement  for  airborne  high  energy  lasers  is  one 
in  which  the  working  medium  flows  radially  outward  from  a  centrally 
located  and  cylindrically  shaped  source.  The  individual  flow  elements  in 
this  arrangement  are  a  bank  of  radial  expansion  nozzles,  an  annular  laser 
cavity,  and  a  supersonic  radial  diffuser.  This  radial  flow  system  is,  in 
particular,  characterized  by  an  axial  length-to-diameter  ( L/D )  ratio  which 
is  many  times  larger  than  that  of  commonly  known  radial  diffusers  such  as 
those  used  in  turbomachines.  This  large  relative  length  introduces 
complications  for  the  design  of  such  systems.  Particular  difficulties 
in  this  respect  are  encountered  with  the  design  of  the  diffuser. 

This  report  describes  the  design  of  a  cold  flow  radial  flow  system 
for  investigating  radial  diffusers  for  laser  applications.  The  radial 
flow  feature  requires  the  design  of  special  flow  expansion  units,  which 
in  the  present  case  consist  of  "ring"  nozzles  (Figure  1).  In  the  first 
part  of  this  report,  a  special  design  procedure  for  these  nozzles  is 
developed.  To  produce  the  radial  flow  field  applicable  to  radial  flow 
lasers  the  expansion  must  take  place  in  a  multiplicity  of  nozzles.  The 
subdivision  of  the  flow  into  small  units  introduces  flow  losses  which 
can  no  longer  be  accounted  for  by  a  simple  expansion  efficiency.  In  the 
second  part  of  the  report,  a  careful  loss  analysis  is  presented  for  an 
analytic  prediction  of  the  expansion  nozzle  performance.  This  prediction 
is  a  necessity  for  the  present  test  rig  design  which  does  not  allow  an 
experimental  performance  determination.  In  the  final  part  of  the  report, 
some  basic  considerations  are  given  to  the  design  of  a  multi-channel 
radial  diffuser. 


AFWAL-TR-80-3028 


SECTION  II 
RING  NOZZLE  DESIGN 

1.  DESCRIPTION  OF  THE  NOZZLE  SYSTEM 

The  radial  outflow  expansion  nozzles  are  made  up  of  ring  shaped 
elements  forming  the  walls  of  the  nozzles.  They  are  referred  to  here 
as  "ring  nozzles".  Figure  1  gives  a  scheme  of  a  ring  nozzle  system. 

The  perforated  tube  in  the  center  of  the  system  is  a  device  to 
provide  uniform  flow  conditions  along  the  axis  of  the  system.  By 
maintaining  a  sufficient  low  porosity  of  the  tube,  the  pressure  in  the 
tube  can  be  kept  high  and  the  velocities  in  the  tube  low.  The  necessity 
for  such  a  supply  tube  depends  primarily  on  the  1 ength-to-diameter  ratio 
of  the  nozzle  system,  i.e.,  the  larger  the  number  of  the  ring  nozzles  the 
higher  the  pressure  must  be.  This  perforated  tube  arrangement  is  a 
specific  feature  of  the  present  model. 

For  an  actual  laser,  either  due  to  heat  addition  in  the  nozzle  system 
(chemical  laser)  or  due  to  very  high  expansion  Mach  numbers  (gas  dynamic 
laser)  noz2le  configurations  result,  which  in  general  do  not  require  the 
use  of  a  high  pressure  supply  tube  to  arrive  at  uniform  flow  conditions. 
The  present  experimental  investigations  deal  with  the  comparatively  low 
Mach  numbers  of  the  chemical  lasers  (M  =  2  to  3),  use  cold  air  as 
operating  medium,  i.e.,  no  heat  addition,  and  contemplate  testing  large 
length  to  diameter  systems  (Jj-  -*■  3).  All  these  conditions  make  the 
perforated  supply  tube  an  absolute  necessity. 

First,  equations  will  be  derived  which  relate  the  nozzle  dimensions 
indicated  in  Figure  1  to  the  fluid-dynamic  and  thermodynamic  conditions 
of  the  flow  system.  These  equations  are  then  used  for  the  preparation 
of  ring  nozzle  design  charts. 
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2.  DERIVATIONS  OF  THE  RING  NOZZLE  DESIGN  EQUATIONS 

The  following  magnitudes  enter  the  ring  nozzle  design  (for  geometric 
data  see  Figure  1 ) : 

Flow  conditions : 

Mach  number  M 
Reynolds  number  Re*  (at  throat) 

Mass  flow  rate  W 
Thermodynamic  conditions: 

Plenum  pressure  Pq 

Plenum  temperature  Tq 

Ratio  of  spec,  heats  y 

Gas  constant  R^ 

Absolute  viscosity  y*  (at  throat) 

With  above-1  isted  magnitudes,  five  equations  can  be  formed  with  R/r 
and  H/h  appearing  explicitly  in  all  equations.  Using  these  two  ratios 
as  coordinates  of  a  graphic  display  of  the  relations,  all  magnitudes 
necessary  for  a  ring  nozzle  design  can  be  represented  in  one  diagram 
(design  chart)  as  shown  in  Figures  2  through  5  for  a  number  of  different 
gases.  The  five  relations  are  derived  below.  They  are  listed  in  terms 
of  the  pertinent  curve  parameter. 

a.  Mach  number 

From  continuity 


Geometric  conditions: 
R 
H 

a 

R/r 

H/h 

R/H 

L/(2R) 


A  =  2RttH  =  R^H 
A*  2r7rh  r*h 


3 
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with  the  well  known  relation  for  the  area  expansion  ratio  in  isentropic 
flow  we  obtain  the  first  equation 


(R)  (H)  =  1 

l  W  VtJ  M 


1  + 


IiIm2 


Y+i 

2 


Y+l 

2T^TT 


(2) 


For  given  gas,  i.e.,  for  a  given  y*  this  relation  can  be  plotted  in 
the  design  chart  with  the  Mach  number  as  parameter. 


b.  (R/H)  tga 

From  Figure  1  we  find 


tga 


H-h 

2 

R-  r 


1  - 


h 

H 


Z(~  -  L) 
^  H  H; 


with  some  rearrangement  we  obtain  the  second  design  equation 


(3) 


H 

h 


R  r\  R  +  R 

—  2  77  tga  —  - 

r  H  y  r 


(4) 


This  relation  is  plotted  in  the  design  chart  with  (R/H) tga  as 
parameter . 


c.  (R/W)Re* 


Continuity  yields 


with 


W  =  A*p*v* 


A*  =  2r  7T-h 


Re* 


v*-h»p* 

v* 


(5) 

(6) 

(7) 
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Equation  5  becomes 


W  =  2rTip*Re* 


(8) 


multiplying  both  sides  by  R  we  obtain  as  third  equation  after  some 
rearrangement 


7  =  2-np* 


[5  -*] 


(9) 


For  a  given  gas,  i.e.,  a  given  p*,  the  expression  (R/W)Re*  is 
proportional  to  R/r  and  can  be  plotted  along  the  R/r  axis  in  the 
design  chart. 


d.  (P0H)/(V^Re*> 


The  continuity  relation  given  by  Equation  5  can  be  developed  in  a 
different  way  by  using  the  equation  of  state 


P* 


T  P* 
o 


V*  RGT0  T*P 


(10) 


and  the  equation  for  the  sonic  speed 


V*  = 


(ii) 


Using  these  two  equations,  Equation  5  becomes 

I  p*  -\  Y 
W  =  A* 1  K 


po  VT*  RG 


(12) 


The  expression  in  brackets  is  a  unique  function  of  the  gas  properties. 
It  will  be  brought  into  a  convenient  form  for  evaluation.  Considering 
sonic  conditions  in  the  throat 


T* 


1  + 


Y-l 


(13) 


s*  ■  <i  ♦ 


(14) 
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the  expression  in  brackets  can  be  written 


With  this  substitution  and  using  Equations  6  and  8,  Equation  12  can  be 
written 


y*Re*  =  h  •  4>  pL 


(16) 


multiplying  both  sides  of  Equation  16  by  H  we  arrive  at  the  fourth  design 
equation 


H 

h 


H_ 

Re* 


4 

b* 


(17) 


The  expression  in  brackets  is  proportional  to  H/h  for  a  given  gas, 
and  can  be  plotted  along  the  H/h  axis  in  the  design  chart.  The  magnitude 
$/i/is  a  unique  function  of  the  working  medium. 


e.  (PExH  r)/(  To  W) 


Using  Equation  6  and  dividing  both  sides  of  Equation  12  by  H  and  P 


this  equation  becomes 


Pc  H 
Ex 


2rTih  P* 
H  P 


*  /V71 

oV 


VT^PEx 


Ex 


(18) 


Using  Equation  15  we  arrive  after  some  simple  transformation  at  the  fifth 
equation 


H 

h 


H  r  P 


Ex' 


w  VTo1 


(19) 
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For  a  given  gas,  the  ratio  P0/PEx  is  a  sole  function  of  the  nozzle  exit 
Mach  number 


(20) 


According  to  the  area  relation  for  a  supersonic  nozzle  (Equation  2), 
Mach  number,  M,  in  Equation  20  and  consequently  the  pressure  ratio  in 
Equation  19  are  both,  for  a  given  gas,  sole  functions  of  the  coordinate 
values  of  the  design  chart.  Thus,  Equation  19  can  be  entered  into  the 
design  chart  with  the  expression  in  brackets  appearing  as  parameter,  0 
being  also  a  unique  function  of  the  gas  properties. 


3.  NUMERICAL  EVALUATION  OF  THE  DESIGN  EQUATIONS 

The  equations  derived  above  have  been  evaluated  for  four  different 
gases  and  the  results  are  plotted  in  Figures  2  through  5.  The  following 
consistent  set  of  dimensions  has  been  used  for  the  evaluation  of  the 
design  equations: 

W  Kg/sec 

H  m 

2 

P  Newton/m 

R  m 

The  required  properties  of  the  here  considered  gases  are  listed  in 
Table  1.  Included  in  the  table  are  expressions  which  occur  in  Equations 
9,  15,  17,  and  19  and  are  derived  from  the  properties.  The  listed  abso¬ 
lute  viscosity  values  apply  only  to  near  ambient  gas  temperatures. 


T  °K 

Rg  Newton  m/(Kg  °K) 

2 

u  Newton  sec/m 


4.  DISCUSSION  OF  THE  RING  NOZZLE  DESIGN  CHARTS 

The  design  charts  in  Figures  2  through  5  apply  to  four  different  gases. 
A  point  in  any  of  the  charts  fixes  the  geometric  design  of  a  ring  nozzle 
unit  for  the  indicated  gas  as  working  medium.  To  fix  the  design  point 
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itself,  either  its  coordinates,  two  curve  parameters,  or  one  of  its 
coordinates  and  one  parameter  must  be  given.  Examples  for  designing 
ring  nozzle  systems  will  be  shown  in  Section  II. 

Two  limiting  conditions  for  the  ring  nozzle  design  can  be  directly 
recognized  from  the  design  charts. 

a.  For  a  multi-ring  nozzle  system,  the  value  of  H/h  should  not  be 
much  lower  than  2;  otherwise  the  subsonic  approach  Mach  number  to  the 
nozzle  becomes  to  high.  This  is  obvious  from  Figure  1  which  shows  that 

in  a  stack  of  nozzles  the  subsonic  flow  cross  section  available  is  limited 
by  the  value  of  H/h.  For  instance,  for  H/h  =  2,  the  subsonic  Mach 
number  at  the  inlet  to  the  nozzle  ring  can  never  be  smaller  than  0.3 
in  case  of  air.  To  obtain  uniform  flow  conditions  in  the  ring  nozzle 
system,  a  low  subsonic  approach  Mach  number  is  desirable  or  extra 
measures  must  be  taken  to  provide  a  uniform  flow  approaching  the  nozzle 
throat  (Section  II). 

b.  Limiting  conditions  are  also  introduced  by  the  values  which  are 
desired  for  the  parameter  (R/H)tga.  Values  for  a  should  be  about  10°, 
i.e.,  tg a  should  be  about  0.2.  R  is  in  general  several  times  larger 
than  H,  say  4  times.  Thus,  ( R/H )  tga  will  be  typically  about  0.8.  Thus 
the  upper  left  hand  and  lower  right  hand  corners  of  the  design  charts 
are  generally  undesirable  regions  for  ring  nozzle  operations. 

In  general,  the  charts  show  that  nozzles  for  low  supersonic  Mach 
numbers,  say  at  least  below  2,  are  difficult  to  design.  Comparing 
Figure  2  to  Figure  5,  where  each  figure  is  made  up  for  a  different  gas, 
one  recognizes  that  the  difficulties  are  eased  with  gases  having  a  low 
y  such  as  CO^  or  Freon. 

In  principle,  there  is  a  ring  nozzle  possible  with  a  =  0°  and  con¬ 
sequently  H/h  =  1,  i.e.,  the  nozzle  walls  are  parallel.  As  the  charts 
show,  very  large  radius  ratios  R/r  are  necessary  to  obtain  high  Mach 
numbers  in  this  case.  Such  a  parallel  wall  expansion  system  still  requires 
a  throat  section  for  the  transition  from  subsonic  to  supersonic  flow. 
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For  a  large  axial  1 ength-to-diameter  expansion  system  where  the  flow 
must  be  subdivided  into  small  expansion  units  to  assure  a  uniform  approach 
flow  to  the  system,  the  parallel  wall  arrangement  affects  the  ring  nozzle 
design  in  a  special  way.  Where  the  throat  contour  changes  to  the  parallel 
wall  contour,  a  subdivision  of  the  flow  is  no  longer  necessary  since 
parallel  walls  in  the  flow  serve  no  useful  purpose  in  this  case.  Thus, 
the  ring  nozzle  system  or  any  other  type  of  subdivision  is  reduced  to  an 
array  of  very  short  nozzles.  In  the  extreme,  the  throat  section  may  be 
only  a  perforated  tube.  Due  to  the  inherent  requirement  for  a  throat  the 
angle  at  as  defined  in  Figure  1  can  never  be  zero  in  a  multi-nozzle 
arrangement . 

The  contour  of  a  nozzle  ring  between  throat  and  exit  is  not  specified 
in  this  design  procedure.  The  angle  a  relates  only  the  location  of  the 
exit  edge  to  that  of  the  throat.  One  can  think  of  a  nozzle  contour  where 
the  flow  expansion  is  made  very  strong  immediately  after  the  throat,  as 
actually  done  for  short,  highly  efficient  nozzles,  with  the  outer  portion 
of  the  nozzle  rings  ending  into  parallel  walls.  Similar  to  the  reduced 
nozzle  system  explained  above,  these  parallel  portions  could  again  be 
eliminated  without  affecting  the  expansion  process.  Thus,  the  design 
charts  give  the  nozzle  Mach  number  at  a  radius  derived  from  the  angle  i 
and  not  necessarily  given  by  the  nozzle  exit  edge. 

5.  DESIGN  EXAMPLES 

For  the  experimental  flow  model  for  which  the  present  design  study 
has  been  initiated  as  mentioned  in  Section  I,  a  Mach  number  2.3  and  a 
Mach  number  3.3  nozzle  system  for  air  as  operating  medium  were  designed. 
The  design  process  followed  in  these  cases  will  be  given  here. 

a .  M  3. 3  nozzl e  system 

Fixed  conditions  for  this  system  were: 

L/(2R)  =  3  (given  condition) 

^  =  9  Kg/sec  (limit  of  air  supply) 

Tq  =  293°K  (test  stand  design) 
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Ml  Ik  N0-3U2K 


if  1 1  .1 1>  W*  *  oihJi  t  ions  were  : 

-  In 

-  u.ul/7  m  (from  supply  tube  design) 

M  1 .  4 

u./  Kg'/ 1  in*"  -  1.962  •  10^  Newton/m^ 

•i  *  n.ijUb  iii  (large  enough  for  flow  survey) 

Re’  hr' 


lx,  a  f  1 1  t  approach  we  try: 

/R  -  1 . 75"  -  0.0444m 

R/r  -  1  .  7b 

,  !*<  H  0.005m 

i  !  Lin.  R/H  4,44 

1 1 ,  !  (R/H;tg,  D./M 

i  f  it  in,  the  number  of  ring  nozzles 

l  -  n  .  h 

) 1  -  L/(2kj  -  3 

2R  •  i 

"  H 

r  -  i  ,?i. light  number'  we  »  house  n 

.St  wt-ijnt  Mow  per  nn/ile  unit  is 

W. 


1.0444  •  1 
0.005 


/  7 


W  to.t  *  ,  0.131  Kg/sec 


.-.-f  g«*  ’  ♦  rn-  i  ne  the  curve  parameter 

,  u  4 

'  f  x  [j.DOb  •  0.0/54  •  0./  ■  9 . M  •  10 

Ia-VT"  0.33  1  •  17.1  •  2 

7  0 

'»  m,i  Millie  and  R/r  1.75  we  read  from  the  chart  in 
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The  actual  dimensions  chosen  for  the  engineering  design  of  the  nozzle 
system  were  the  following; 


2R  =  1 .75" 

-  44.45mm 

-  .04445m 

2r  =  1 .00“ 

=  25.40mm 

=  .0254m 

H  =  .192" 

=  4.877mm 

=  .004877m 

h  -  .059" 

-  1 .499mm 

=  .001499m 

.  =  10" 

Only  H  is  slightly  different  from  the  value  used  in  the  first  approach 
From  these  chosen  values  follows: 


R/r  - 

1  .75 

A/A* 

=  5.695 

H/h  - 

3.254 

R/H 

-  4.557 

M 

3.  31 

(R/H) tg  . 

-  0.808 

The  other  magnitudes  occurring  in  the  design  chart  become 


p 

Ex  •  r  H 


0.?? 


"  Re*  --  1.547  104 

W 


P  H 
0 

yF  Re* 

V  o 


1.45  ■  103 


Newton  sec 
Kg  "K 

m  sec 

*g 


Newton 
m  °K 


Derived  values  of  interest  are: 


=  0.22 


60701 

m  Kg 


=  0.6188  -s-e-c 
in  Kg 


/ 
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(  2 

tor  example,  for  =  0.2  Kg'/m 
it  is  W  =  oTifs  =  0,323  K9/sec 
For  the  nozzle  throat  Reynold  number,  we  obtain 
Re*  -  1.547  •  1 04  •  5 

K 

=  2.25  •  105 

The  nozzle  system  in  Figure  1  is  drawn  in  proportion  to  the  designs 
presented  above. 

h.  M  2 . 3  Nozzl e  System 

This  system  was  intended  to  have  the  same  outside  diameter  and 
nozzle  spacing  as  the  high  Mach  number  nozzle  system. 

Thus,  the  following  magnitudes  are  fixed: 

2R  =  0.0444  m 
H  =  0.00488  m 
R/H  -  4.557 

Given  by  the  test  facility:  T  =  293°  K 

7  o 

For'  magnitudes  which  can  be  freely  chosen  the  following  desirable  ranges 
were  es tabl i shed : 

t  =  10°  to  15° 

R/r  •  1  .2 
H/h  >  1.5 
M  2.3 

It  follows  that  (R/H)  tg  i  =  0.803  to  1.222 
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Those  ranges  fix  an  area  in  the  design  chart  within  which  the  design 
point  may  be  selected,  as  shown  in  Figure  3.  The  following  coordinates 
were  chosen  for  the  point: 

R/r  =  1.25  H/h  =  1.7 


The  values  of  the  curve  parameters  for  the  chosen  point  are: 


(R/H)tga  =1.03 
M  =  2.265 


P 


M 


0.566 


Newton  sec 
Kg  °K 


Other  derived  values  are: 

tg.  =  0.2312 
.  =  13.02° 

-£■  =  1  56263  S6C  =  1.593 

m  Kg  crn  Kg 

2 

f°r  =  0.5  Kg'/cm  it  is  W  =  0.314  Kg/sec 

A  low  Mach  number  ring  nozzle  system  designed  according  to  conditions 
above  is  shown  in  Figure  6.  As  indicated  before,  due  to  the  low  H/h 
value  of  1.7,  the  subsonic  approach  Mach  number  is  near  .4  and  even  higher 
at  the  smaller  radius  where  the  screen  is  positioned.  This  screen  serves 
as  flow  equalizer  for  the  flow  approaching  the  actual  expansion  nozzle. 
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SECTION  III 

RING  NOZZLE  PERFORMANCE  ANALYSIS 

1.  GENERAL  CONSIDERATIONS 

How  expansion  in  a  supersonic  nozzle  is  normally  a  very  efficient 
proi.es is .  For  instance,  for  a  Mach  number  3,  axis-symmetric  nozzle 
etti.  iencies  near  99  are  possible.  However,  for  the  radial  outflow 
./stem  considered  here,  lower  efficiencies  must  be  expected  due  to  certain 
geometric  constraints  imposed  on  the  system.  Laser  applications  require 
that  the  axial  length  of  the  nozzle  system  is  several  times  larger  than 
it*,  diameter.  Such  shape  can  only  be  obtained  by  subdividing  the 
r  a  pa  ns  ion  '/stent  into  many  small  units.  In  the  present  case,  this  is 
done  Ly  i  ,eries  of  nozzle  rings  (Figure  1).  This  subdivision  reduces 
’mi  i>.  ter  ist  ic  length  dimensions  and  consequently  the  characteristic 
-•/h.;  l  number  of  the  system.  Thus,  for  the  same  flow  rate,  the  wall 
.  r.  i  <  >  n  coefficient  for  the  radial  outflow  system  must  be  higher  than 
•  >u  *  *01  an  a x i s  - syinme tr i c  nozzle. 

At  the  exit  of  t he  rnulti  unit  expansion  system,  the  flow  experiences 
j  -ji  Men  flow  t  ross- sect  ion  enlargement  due  to  the  finite  thickness  of 
♦fn:  nozzle  exit  edges.  This  is  another  cause  for  flow  losses.  Finally. 
*’(  t’uw  in  a  radial  outflow  system  has  an  inherent  divergence  which 
,  *  f  j  *  i  «  t;i  ri  momentum  loss,  if  the  flow  has  to  be  turned  into  the  discrete 
•Hr*,  tines  uf  the  channels  of  the  radial  diffuser. 

I  he  expansion  losses  can  be  determined  fairly  closely  by  analysis. 

Ihis  is  important  for  the  evaluation  of  experimental  results  with  radial 
,?f!ow  s/stems.  The  end  conditions  of  the  expansion,  which  are  the 
n  !* rondiMnn  tor  the  diffuser  cannot  clearly  be  determined  in  the 
•  imcnt'. ,  ma  i  n  1  /  due  to  the  presence  of  the  nozzle  wakes  in  the  flow. 

Mu  pace  between  the  nozzle  exits  and  the  diffuser  inlet  (laser 
-,/}♦,/  is  not  readily  accessible  to  i  nstrumenta  tion  for  a  complete  flow 
j  r  v1  *  ■  / .  Analytic  predictions  which  account  for  possible  losses  during 
t-,\  in,  ion  are  therefore  more  reliable  than  limited  flow  measurements. 
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In  the  following,  an  analytic  process  for  calculating  the  performance 
of  a  ring  nozzle  system  accounting  for  the  expansion  losses  indicated 
above  will  be  derived.  The  analysis  relates  the  flow  conditions  between 
significant  flow  cross  section  to  each  other  with  the  conservation  laws 
for  mass,  momentum,  and  energy.  The  resulting  equations  are  given  in  the 
form  of  force  equations  which  allow  a  ready  addition  of  the  wall  friction 
forces  in  terms  of  the  common  pipe  friction  coefficient.  This  method 
has  been  used  previously  for  determining  ejector  performances  (Refer¬ 
ences  1,2).  Equations  developed  for  the  ejector  case  can  be  used  here 
directly  by  setting  the  primary  and  secondary  velocity  equal  in  these 
equations.  In  Reference  2,  this  procedure  has  been  already  applied  to 
the  diffuser  case.  The  derivations  provided  below  follow  very  closely 
the  outline  given  for  the  ejector  case  in  Reference  1. 

2.  DERIVATION  OF  THE  PERFORMANCE  EQUATIONS 

A.  Momentum  Equation 

The  momentum  equation  applied  to  the  adiabatic  flow  in  a  duct,  where 
inlet  (1)  and  exit  (2)  may  have  different  cross  sectional  areas,  can  be 
written  in  the  following  form: 


m-vl  =  m‘v2  +  FW  +  Fp  +  Ff 


(21) 


Fw  in  this  equation  is  the  wall  pressure  force  acting  between  the  wall 
and  the  flow  in  flow  direction  in  a  non-constant  area  duct 

Ai 

Fw  J  (p  1  -  pw}  dA  (22) 

A2 


where  Pu  is  the  local  wall  pressure  and  (dA)  a  differential  wall  area 
w  w 

projected  in  flow  direction. 
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F 

P 

central 


is  the  force  acting  in  the  fluid  in  flow  direction  within  the 
region  where  no  wall  pressure  forces  act 


FP  ■  Mpi 


(23) 


F^  in  Equation  21  is  the  wall  friction  force  acting  on  the  flow  along  the 
duct  wall.  This  force  may  be  expressed  in  the  form  used  for  normal  pipe 
friction 


F 


f 


L 

D 


(24) 


where  a  judicious  choice  must  be  made  for  the  reference  crossection  A. 

The  friction  coefficient,  c^.,  in  this  equation  applies  to  fully  developed 
flow  (Section  IV).  The  fluid  pressure  forces  and  the  friction  forces 
will  be  discussed  below  in  more  details. 


(1)  Fluid  Pressure  Forces  FtJ  and  F 

W  p 

Both  forces  have  the  same  cause,  which  is  the  pressure  change  in  flow 
direction.  They  are  best  treated  together.  Following  procedures  given 
in  Reference  1,  Equation  22  is  written  in  a  form  which  facilitates  the 
integration  of  this  equation  by  introducing  a  shape  factor  i.  All 
equations  are  also  written  in  a  form  consistent  with  the  ejector  case 
For  this  reason  Equation  23  does  not  appear  explicitly  in  che  following 
derivations.  In  this  way,  the  wall  pressure  forces  become 

Pi  ~  P  9 

Fw  =  (A1  -  A2)  x—r  --  i  (25) 

Then  for  the  total  fluid  pressure  forces  one  can  write 

rwtp  •-[(/>,-  A2)i4^t  *A2  (p,  -  P2)J  (26) 

by  using  the  abbreviation: 

a2 

t  =  r  (27) 

Ai 
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Equation  26  can  be  transformed  into 

FW+P  =  ^  p2  '  P1 ^A2^2t  +  1  "  2^ 


1 


+  1-4 


Wri ti ng 

1  =  2 1  '  '  "  2 
the  final  form  of  the  fluid  pressure  forces  becomes 

W1  •  1  ■  »  ■  Vff- 11 


(28) 


(29) 


(30) 


(2)  Wall  Friction  Force 

The  use  of  the  pipe  friction  coefficient  for  expressing  friction 
forces  requires  the  selection  of  a  reference  cross  section  together  with 
the  flow  velocity  in  this  cross  section.  For  a  normal  pipe,  the  reference 
conditions  are  those  at  the  inlet  of  the  pipe.  As  will  be  shown  in 
Section  IV,  in  the  present  case  the  conditions  at  the  duct  exit  are  a 
particularly  useful  choice. 


Considering 


2  2 
p2  v2  =  P2  Y  ^2 


(31) 


we  can  write  for  the  friction  forces  with  Equation  24 


,2  L 


Ff  =  A2  ‘  Cf  '  P2  ‘  Y  ‘  M2  2D 


(32) 


(3)  Final  Momentum  Equation 
Expressing  the  flow  momentum  in  Equation  21  by 

A  m2 

m  *  v  =  A  •  p  *  y  •  M 


(33) 
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and  using  above  developed  expressions  for  fluid  pressure  forces  and 
wall  friction  forces.  Equation  21  becomes 

A1P1yM^  =  VP2-Y-M2  +  T-t-PrAl(^  -  D  +  A2cfP2vM2  ^  (34) 

This  equation  can  be  readily  transformed  info 


+  T 


M2'2ET  +  1  ^ 


+  T 


B 


(35) 


Equation  35  is  valid  for  any  adiabatic  change  of  state  in  a  duct  with 
varying  cross-section  areas  such  as  flow  expansion  in  a  nozzle,  sudden 
change  in  flow  cross-section,  or  supersonic  diffusion  in  a  constant  area 
duct.  The  differentiation  between  these  various  possible  changes 
appears  only  in  the  wall  force  parameter  x.  The  determination  of  t  will 
be  shown  in  Section  III. 


Equation  35  yields  the  pressure  ratio  across  the  duct 


P2 

P1 


B 


2  CfL 

ym2(^+1) 


+  T 


(36) 


The  duct  exit  Mach  number  on  the  right  side  is  still  unknown.  We 
obtain  this  Mach  number  from  the  requirement  for  conservation  of  mass 
and  energy  in  the  flow,  as  shown  below. 

b.  Mass  and  Energy  Equation 
From  continuity,  follows 

m  =  Vj  *  A^  •  =  v^  •  A^  *  ?2  (37) 
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wi  th 


a 


=  >/  V  •  R  •  T 


we  can  write 


or 


P_ 

RT 


"i .PHo 

*  P2\r0  T1 


For  an  adiabatic  change  of  state,  it  is 


2 

2 


Then,  Equation  42  can  be  written 


m20 


+  M*) 


1/2 


With  P-|/P2  from  Equation  36,  we  arrive  at 


m20  +  n\) 

2  ^  f 

y  +  1 )  +  T 


1/2 


'38) 

(39) 

(40) 


(41  ) 


(42) 


(43) 


(44) 


(45) 


The  right  side  of  this  equation  is  made  up  of  magnitudes  strictly 
given  by  the  inlet  conditions  and  the  duct  geometry  in  terms  of  t. 
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The  left  side  contains  only  magnitudes  given  by  the  exit  conditions  and 
the  factor  .  which  is  a  constant  for  a  particular  kind  of  duct  flow. 
Equation  45  can  be  solved  for 


C  my  -  Y  +  1 


(46) 


where 


1  + 


CfL 

2D 


(47) 


«  =  2 


YE2c 


(47a) 


The  positive  root  in  this  equation  gives  the  subsonic  solution  for  the 
exit  Mach  number  and  the  negative  root  the  supersonic  solution  if 
compatible  with  the  existing  conditions. 


3.  WALL  FORCE  PARAMETER 

a.  Expansion  Nozzle 

For  zero  wall  friction,  Equation  46  applied  to  an  expansion  nozzle 
must  yield  the  same  result  as  the  common  thermodynamic  relations  for  an 
isentropic  expansion.  From  this  requirement,  the  factor  r  can  be 
determined  for  the  case  of  an  expansion  nozzle.  With  a  T-value  found  in 
this  way  for  a  particular  case  of  expansion  nozzle,  Equation  46  allows 
one  to  add  friction  forces  to  this  expansion  process  and  find  the  new 
exit  conditions.  This  procedure  to  calculate  flow  expansion  with  friction 
is  an  approximate  method  since  wall  friction  may  change  the  pressure 
distribution  along  the  expansion  nozzle  and  thus  the  value  :  of  the 
process.  There  are  indications  that  this  influence  on  t  is  in  practical 
cases  small  enough  in  comparison  to  the  direct  influence  of  friction  on 
the  flow  that  the  change  of  a  can  be  neglected.  A  possible  approach  for  elim¬ 
inating  this  influence  is  to  calculate  the  flow  expansion  in  segments  small 
enough  for  the  pressure  distribution  to  become  a  straight  line  independent 
of  the  presence  of  wall  friction.  The  factor  x  becomes  a  constant  in  this  case. 
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To  determine  i  for  an  isentropic  expansion  process,  we  solve 
Equation  34  for  r  with  wall  friction  zero 

2  2 
P1yM1A1  -  p2yM2A2 

(p2  -  Pj)  A^ 


(48) 


or  somewhat  rearranged 


2  A1 

m:  — 

1  A?  p2 


(49) 


The  area  ratio  and  the  pressure  ratio  in  this  equation  are  determined 
from  the  well  known  relations  for  isentropic  expansion 


A 

A 


1 

2 


1  + 

2  nl 

1  +TM2J 


•y+1 

2CHT 


(50) 


P2 


1 


1 


v-1  u2' 
+  2  M1 


(51) 


For  a  given  gas,  a  given  inlet  Mach  number  and  a  given  nozzle  area 
ratio,  t  can  be  determined  from  Equation  49  with  the  help  of  Equations  50 
and  51.  The  value  of  t  changes  strongly  with  the  expansion  conditions 
and  must  be  determined  for  each  individual  case. 

b.  Sudden  Crossection  Increase 

For  a  sudden  cross  section  increase,  the  wall  forces  are  given  by: 

FW  =  (A2  ■  A1 )  (p2  ■  P1 ^  (52) 


l 
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Comparing  this  equation  with  Equation  25  one  finds  that  Equation  25  is 
identical  with  the  above  equation  for  an  i-value  of  2.  With  this  value, 
we  obtain  from  Equation  29: 


T  =  \  (53) 

c.  Constant  Area  Duct  (Supersonic  Diffuser) 

For  a  constant  area  duct  with  either  subsonic  or  supersonic  flow, 
the  factor  i  must  be  zero,  since  no  wall  forces  in  flow  direction  can 
appear. 

With  Equation  29 


t  =  1 


(53a) 


4.  FLOW  DIVERGENCE 

In  the  present  radial  diffuser  system  (Figure  7),  the  radial  flow 
must  be  redirected  into  individual  diffuser  channels  as  schematically 
shown  in  Figure  8.  Since  this  redirection  of  the  flow  must  be  accom¬ 
plished  in  a  more  or  less  abrupt  way,  it  is  a  source  for  flow  losses. 
These  losses  can  be  accounted  for  by  applying  the  law  of  conservation 
of  momentum  to  the  flow.  For  the  ratio  of  the  momentum  of  the  flow  in 
direction  of  the  diffuser  channel  to  the  total  momentum  of  the  radial 
flow  portion  entering  the  diffuser  channel,  we  can  write  (for  notations 
see  Figure  8). 

/x# 

.  .P-Y-M2chancos2e  dx 
1  /  p-y.M  d(r-8) 

For  small  angles,  6  (<10°)  the  following  approximations  apply: 

X  =  r  •  8 

3_  =  >L 

P'O  xo 

r2 

cose  ~  i  - 
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Then  the  integration  yields: 

I  «2  ,4 

chan  ‘  o  *  o 

I  1  ’  3  20 

g2 

•*  cos3_  +  o 
0  6 

For  the  diffuser  entrance  Mach  number,  we  can  then  write: 


Resulting  losses  are  in  general  very  minor  for  the  cases  considered  here. 
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SECTION  IV 

WALL  FRICTION  LOSSES  IN  AN  EXPANSION  N077LE 

1 .  GENERAL  CONSIDERATIONS 

In  applying  the  common  concept  of  a  pipe  friction  coefficient  to 
a  duct  of  non-cyl indrical  shape  such  as  an  expansion  nozzle,  two  obvious 
problems  arise.  One  is  the  choice  of  the  proper  flow  condition  to  be 
entered  in  Equation  24;  the  other  one  is  the  choice  of  the  proper  duct 
1 ength-to-diameter  ratio  also  occurring  in  this  equation.  In  the 
following  it  will  be  shown  that  these  problems  have  a  satisfactory 
solution  and  that  there  is  no  formal  difficulty  determining  the  wall 
friction  forces  for  an  expansion  nozzle  as  long  as  the  friction  coefficient 
is  a  given  constant. 

2.  CHARACTERISTIC  FLOW  CONDITIONS  FOR  WALL  FRICTION  LOSSES 
Introducing  in  Equation  32  the  following  substitution 

°2  * 

and  by  introducing  pQ  and  the  square  root  of  A*  into  the  equation  we 
obtain: 


For  a  given  gas,  the  dimensionless  expression  in  brackets  is  a  unique 
function  of  the  Mach  number  prevailing  in  the  exit  crossection  A^,  in 
accordance  with  the  relations  for  isentropic  expansion: 
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This  approach  is  not  quite  accurate,  since  in  the  presence  of  wall 
friction  the  expansion  is  no  longer  isentropic  and  Equations  56  and  5/ 
need  a  correction.  This  could  be  accomplished  by  introducing  an  estimated 
polytropic  expansion  efficiency.  However,  compared  to  the  inaccuracies 
involved  in  determining  the  friction  coefficient  itself,  this  effect  is 
small  and  neglected  here. 

The  magnitudes  in  Equation  55  outside  the  brackets  are  constants  for 
a  given  expansion  nozzle,  except  the  length  L  which  can  be  chosen  freely 
If  we  make  L  very  short  in  comparison  to  the  total  length  of  the 
expansion  nozzle,  we  obtain  the  friction  force  acting  on  a  short  and 
nearly  cylindrical  segment  of  the  expansion  nozzle  at  the  exit  cross- 
section  Ap .  By  applying  this  approach  to  other  cross  sections  along 
the  expansion  nozzle,  we  obtain  the  distribution  of  the  local  friction 
force  along  the  nozzle.  Figure  9  gives  the  distribution  of  the  local 
friction  force  with  the  constant  magnitudes  in  Equation  55  set  unity,  as 
a  function  of  the  square  root  of  A^/A*.  Plotting  against  the  square 
root  of  the  area  ratio  makes  the  distribution  curve  identical  to  the 
one  applicable  to  a  round  expansion  nozzle  with  straight  contours  as 
drawn  in  the  plot.  On  the  abscissa,  the  Mach  number  applicable  to  the 
area  ratio  is  also  entered  in  Figure  9. 

Within  the  accuracy  of  the  chosen  friction  coefficient,  the  area 
underneath  the  distribution  curve  from  the  inlet  to  the  exit  of  the 
expansion  nozzle  is  a  direct  measure  for  the  complete  friction  force 
acting  on  the  nozzle  wall.  With  the  help  of  this  integrated  value,  the 
validity  of  the  simplified  approach  taken  in  the  momentum  equation  in 
Section  III  can  be  checked.  In  Equation  32  the  reference  flow  conditions 
are  those  at  the  exit  of  the  nozzle.  The  L/D  in  this  equation  is 
obtained  by  dividing  the  total  length  of  the  nozzle  by  the  exit  diameter, 
i.e.,  the  implied  assumption  is  that  the  friction  force  is  constant  along 
the  nozzle.  In  Figure  9,  the  distribution  for  this  simplified  approach 
is  then  represented  by  a  straight  line  drawn  parallel  to  the  abscissa 
through  the  exit  point  on  the  true  distribution  curve.  Taking  as  an 
example  the  3.3  Mach  number  nozzle  of  the  present  design  and  assuming 
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an  inlet  Mach  number  of  0.25,  the  areas  beneath  the  straight  line  and  the 
true  curve  are  approxima tely  equal;  i.e.,  both  approaches  give  the  same 
result  in  this  case.  This  coincidence  is  accidental.  However,  the 
distribution  curve  in  Figure  9  is  universally  valid  and  a  correction 
factor  can  be  established  for  any  combination  of  inlet  and  exit  Mach 
number.  These  correction  factors  are  universally  valid  for  any  cross- 
sectional  shape  of  the  expansion  nozzle  (Section  IV. 3). 

With  the  availability  of  a  universally  valid  correction  system,  no 
formal  difficulties  exist  to  obtain  correct  friction  forces  with  the 
simplified  method  introduced  here. 

The  real  difficulties  for  determining  friction  forces  in  supersonic 
expansion  nozzles  lie  with  the  friction  coefficient  itself.  Wall  friction 
coefficients  for  supersonic  flow  are  not  well  established.  Besides  being 
a  function  of  the  Reynolds  number,  they  depend  also  on  the  Mach  number 
itself.  In  addition,  the  flow  in  an  expansion  nozzle  is  never  the  fully 
developed  type  of  flow  normally  assumed  for  standarized  wall  friction 
determi nations . 

The  advantage  of  the  present  method  of  wall  friction  determination 
in  expansion  nozzles  is  that  it  allows  comparative  studies  of  expansion 
nozzle  configurations  and  provides  at  least  approximate  results  in  the 
absence  of  any  other  simple  method. 

3.  NON-CIRCULAR  EXPANSION  NOZZLE 

The  common  relation  for  pipe  friction  as  used  here  assumes  a  circular 
cross  section  for  the  duct  for  which  the  friction  forces  should  be 
determined.  This  restriction  is  only  a  formal  one.  In  the  basic 
derivation  of  the  pipe  friction  relation,  this  restriction  does  not 
even  enter  since  it  deals  with  the  concept  of  the  hydraulic  radius 
which  is  by  definition  independent  of  the  shape  of  the  cross  section. 

The  common  pipe  friction  relation  is  only  a  special  case  of  the  general 
relation. 
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The  basic  relation  for  friction  on  a  wall  is: 

Ff  =  CW  2“  AW 


For  a  cylindrical  duct  of  any  cross  section  shape,  the  wall  surface  is 


*w  ■  U'L 


(59) 


By  writing 


A 


W 


■L-A 


cr 


cr 


(60) 


the  duct  cross  section  area  can  be  introduced  in  the  relation  for  the 
friction  force  in  a  duct.  The  ratio  A^/U  is  the  well  known  hydraulic 
radius  of  the  duct: 


rhyd 


(61) 


With  this  relation  Equation  58  can  be  written: 


F 


f 


(62) 


This  is  the  basic  relation  for  the  friction  force  in  a  duct  of 
arbitrary  cross  section  shape.  Applying  this  relation  to  a  circular 
cross  section  with  diameter  6,  by  writing 


=  ifn  *  £ 

rhyd  46tt  4 


(63) 


Equation  62  becomes 


p*v 


r  =  4  •  c  _ 

f  W  2 
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(64) 


It  is  a  convention  to  set 

4  c  = 

w 

calling  the  pipe  friction  coefficient. 


(65) 
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T  .  ot  Mii'.  |)>j>v  friction  coefficient,  Equation  6?  becomes 

r  !C  ■/.  l.  .  A 

f  f  2  4rhyd  Cr 


(66) 


runs  in  any  relation  previously  derived,  the  friction  parameter: 


•  irL 

C  '  2D~ 


(67) 


in  he  replaced  hy  the  more  general  expression  with  writing  2r^  ^  =  ^^yd 


C-.  Cf  L 


4d 


hyd 


(68) 


Both  equations  give  the  same  result  for  a  circular  cross  section. 

For  a  rectangular  cross  section  which  is  of  particular  concern  for 
a  ring  nozzle  system,  it  is  (with  small  side  a  and  long  side  b) 


nyd  (g  +  1 ) 


(69) 


Equation  66  can  then  be  written: 


F.  =  c,  —  (  r-  +  1 )  A 
f  f  2  a  b  cr 


(70) 


Tilt-  friction  parameter  is  then: 


c'-5r  <!  +  1> 


(71) 


4.  WALL  FRICTION  COEFFICIENT  FOR  AN  EXPANSION  NOZZLE 

The  prediction  of  a  wall  friction  coefficient  for  an  expansion  nozzle 
is  a  complex  process.  For  a  straight  pipe,  the  wall  friction  is  already 
a  function  of  Reynolds  number,  Mach  number,  wall  surface  conditions, 
and  flow  history.  In  case  of  an  expansion  nozzle,  all  these  influence 
parameters  change  to  an  appreciable  degree  from  inlet  to  exit  of  the  nozzle. 
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The  influence  of  the  Mach  number  on  the  friction  coefficient  is  not 
very  well  established.  Since  this  influence  can  be  considered  minor  in 
comparison  to  those  of  the  other  parameters  for  the  accuracies  desired 
here,  it  is  neglected  in  the  present  calculations.  The  flow  history  is  an 
important  factor.  Since  the  flow  approaching  the  ring  nozzles  is  pre¬ 
expanded  in  the  perforations  of  the  supply  tube  (Figure  1)  and  abruptly 
dispersed  before  entering  the  ring  nozzles,  it  is  extremely  turbulent 
and  the  flow  profile  development  is  likely  to  occur  very  fast.  The  large 
pressure  gradients  along  the  flow  in  the  nozzle  also  promote  a  fast 
development  of  the  velocity  profiles  at  each  station  in  the  nozzle.  For 
the  present  calculations,  the  flow  in  the  nozzle  is  therefore  assumed 
to  be  fully  developed  in  a  sense  that  at  any  local  station  along  the  flow 
path  in  the  nozzle,  the  velocity  profile  has  the  character  of  a  fully 
developed  flow. 


The  relative  wall  roughness,  given  by  the  ratio  of  passage  diameter 
over  wall  roughness  height,  was  estimated  to  be  500  on  the  average  for  the 
ring  nozzles  used  in  the  present  investigations.  As  we  will  see  below, 
this  degree  of  roughness  makes  the  Reynolds  number  influence  nearly 
constant  along  the  path  of  expansion  in  the  ring  nozzles.  To  show  this, 
the  change  of  the  Reynolds  number  along  the  path  of  expansion  will  be 
derived  in  the  following. 


One  can  write  for  the  Reynolds  number  using  the  equation  of  state: 


P 


P  •  Y 

VT'T 


(7?) 


or: 


Re.MP  — - 

d  '  U 


This  relation  can  be  extended  to  read: 


(73) 


(74) 
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The  pressure  and  temperature  ratio  in  this  relation  can  be  expressed 
in  terms  of  the  occurring  Mach  number: 


assuming  u  to  be  a  constant* 
With  the  abbrevia tion : 


1 

u 


Equation  75  becomes 


(75) 


(76) 


Re  _ 
d 


•X 


(77) 


For  air,  this  relation  gives  the  following  values: 


M  = 


1  2  3  3.3 

2.25  1.33  0.53  0.40 


Taking  d  in  Equation  77  as  the  local  hydraulic  diameter  of  the 
expansion  nozzle: 


exit  area  ,  _  (2R~)  h  _  h 

circumfer.  of  exit  h  (2R^)  2  2 


(78) 


we  find  for  the  ring  nozzles  used  in  the  present  experiments: 

throat  nozzle  exit 


h 

2 

Re  = 


2.25  0.40 

0.75  •  10'3  m  2.3  •  10‘3  m 

2  -1 05  1.1  •  105 
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From  the  well  known  plot  of  the  pipe  friction  coefficient  versus 
the  Reynolds  number  we  find  a  value  of: 

c  f  =  0.02 

for  the  pipe  friction  coefficient  at  the  throat  as  well  as  at  the  nozzle 
exit.  This  equality  is  caused  by  the  wall  roughness  encountered  here 
which  makes  the  pipe  friction  coefficient  nearly  independent  of  the 
Reynolds  number  in  the  pertinent  region. 
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SECTION  V 

THE  THEORETICAL  PERFORMANCE  OF  A 
RADIAL  OUTFLOW  SYSTEM 

1.  GENERAL  ASPECTS 

In  determining  the  theoretical  performance  of  a  radial  outflow 
system  as  schema ti cal  1 y  shown  in  Figure  7,  only  the  expansion  process 
is  accessible  to  a  realistic  treatment.  The  diffusion  process  is  much 
too  complex  to  allow  any  closer  performance  prediction  than  that  applicable 
to  normal  shock  pressure  recovery.  As  explained  in  Section  III,  the 
present  arrangement  for  investigating  radial  flow  diffusers  relies  on 
the  predi ctabi 1 i ty  of  the  expansion  performance.  In  this  Section,  the 
results  for  the  Mach  No.  3.3  ring  nozzle  expansion  process  will  be 
presented  with  consideration  given  to  the  likely  accuracy  of  the 
calculations. 

For  the  diffuser  portion  of  the  flow  system  only  the  normal  shock  re¬ 
covery  for  the  diffuser  inlet  flow,  which  is  the  same  as  the  expansion 
exit  flow,  has  been  derived.  This  normal  shock  performance  serves  as 
standard  for  evaluating  diffuser  experiments.  The  present  calculations 
allow  the  inclusion  of  wall  friction  losses  in  determining  the  normal 
shock  performance.  Results  given  in  this  Section  demonstrate  the  com¬ 
paratively  small  performance  penalty  to  be  paid  for  subdividing  the 
radial  diffuser  into  many  small  units  exposing  a  large  friction  area. 

An  important  subject  discussed  in  this  Section  is  the  independence 
of  the  present  performance  calculations  from  the  non-uni formi ti es  of 
the  flow.  This  independence  strongly  benefits  the  reliability  of  the 
calculations. 

2.  PERFORMANCE  CALCULATIONS 

Table  2  provides  the  design  conditions  and  the  analytic  performance 
of  the  expansion  and  diffusion  process  for  the  Mach  number  3.3  ring 
nozzle  system  treated  in  Section  II.  The  first  column  lists  the 
individual  steps  in  which  the  calculations  proceed.  The  group  of  the 
next  seven  columns  gives  the  necessary  design  data  together  with 
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immediate  derivatives  such  as  the  wall  force  parameter  or  the  hydraulic 
diameter  of  the  exit  cross  sections.  The  applicable  equations  for 
determining  i  for  each  step  of  the  flow  process  are  listed  in  the  column 
heading.  The  next  three  columns  give  the  results  of  the  performance 
calculations  in  terms  of  expansion  Mach  numbers,  normal  shock  recoveries 
with  and  without  wall  friction,  and  expansion  pressure  ratios.  The 
equations  used  to  obtain  these  performance  data  are  indicated  again  in 
the  column  headings.  For  comparison,  the  performance  results  for  an 
isentropic  expansion  process  are  also  listed  (see  discussion  on  accuracy). 

The  essential  result  for  the  evaluation  of  the  diffuser  experiments 
is  the  normal  shock  pressure  recovery  ratio  14.24  listed  in  column  10 
for  the  diffuser  without  wall  friction.  Experimental  diffuser  pressure 
recovery  ratios  obtained  with  the  present  ring  nozzle  and  cavity 
arrangement  are  to  be  directly  compared  with  above  figure  for  judging 
the  diffuser  performance.  For  supersonic  diffusers  with  a  more  elaborate 
geometry  than  a  constant  area  duct,  this  comparison  figure  is  not 
necessarily  an  upper  limit  of  performance.  In  general,  however, 
obtaining  normal  shock  pressure  recovery  with  any  conventional  type  of 
supersonic  diffuser  is  considered  a  very  good  performance  (Reference  4). 

The  performance  difference  in  Table  2  for  the  diffuser  with  and 
without  wall  friction  is  of  basic  interest  for  applying  the  multi-channel 
diffuser  concept  in  the  present  investigations.  This  concept  facilitates 
the  fl ui d-dynamic  design  of  radial  diffusers  as  explained  in  the  next 
Section.  The  performance  with  wall  friction  that  applies  to  the  condition 
of  a  multi-channel  arrangement  with  a  commonly  required  channel  length 
to  hydraulic  diameter  ratio  of  12,  is  less  than  3-  lower  in  terms  of 
pressure  recovery  as  compared  to  the  frictionless  diffuser,  as  a  close 
representation  of  a  diffuser  with  no  subdividing  walls.  For  completeness, 
it  must  be  mentioned  that  actual  radial  flow  laser  systems  would  also 
incorporate  a  subsonic  diffuser.  Since  subsonic  diffusers  are  not 
considered  a  specific  problem  of  the  present  radial  flow  system,  they 
have  been  omitted  from  these  investigations. 
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i.  ACCURACY  OF  THE  PERFORMANCE  CALCULATIONS 

The  analysis  provides  a  standard  of  comparison  for  experimental 
diffuser  results.  As  indicated  above,  this  standard  is  a  sole  function 
of  the  expansion  exit  conditions.  Two  factors  influence  the  accuracy 
with  which  these  conditions  can  be  predicted.  One  is  the  inherent  non- 
uni  formity  of  the  expansion  flow;  the  other  one  is  the  proper  account  of 
the  flow  losses. 

Due  to  the  previously  discussed  necessity  of  subdividing  the  flow 
expansion  system  into  small  units,  the  expansion  flow  contains  wakes 
which  cause  a  considerable  non-uni formi ty  of  the  flow  entering  the 
diffuser.  The  nature  of  the  supersonic  diffusion  process  allows  a  ready 
account  of  this  non-uni formi ty .  Since  the  supersonic  diffusion  process 
i.  based  on  the  preservation  of  the  flow  momentum,  the  diffuser  inlet 
onditions  need  only  be  interpreted  in  terms  of  the  flow  momentum  entering 
the  diffuser.  For  an  experimental  determination  of  this  momentum,  the 
impact  forces  measured  with  a  pitot  probe  must  be  integrated  over  the 
i'-let  area.  For  the  analytic  determination  which  is  the  only  one  possible 
as  explained  before,  the  flow  momentum  is  an  inherent  result  of  the 
present  tabulations,  which  describe  the  expansion  process  in  terms  of 
the  prevailing  flow  forces  (Equation  21).  The  effective  Mach  number  of 
tne  expansion  process  resulting  from  the  calculations  stands  as  a  direct 
criterion  for  the  integrated  exit  flow  momentum. 

For  the  proper  account  of  the  flow  losses,  it  must  be  considered  that 
the  total  losses  involved  are  approximately  1 0C/  of  the  available  expansion 
energy  (compare  column  10  and  14  in  Table  2).  Two  percentage  points  of 
this  figure,  as  can  be  derived  from  the  numbers  given  in  column  9,  are 
due  to  a  sudden  change  in  flow  cross  section  and  direction  of  the  flow 
;nth  (Section  III).  This  influence  is  a  function  of  the  flow  geometry 
.in si  can  be  readily  accounted  for.  As  a  result,  about  8'  of  the  initial 
expansion  energy  is  lost  due  to  wall  friction  which  can  be  determined 
only  with  a  limited  accuracy.  The  uncertainty  is  mainly  associated 
with  r.h*'  friction  coefficient  to  which  the  losses  are  directly  proportional. 
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An  estimate  of  the  accuracy  of  this  coefficient  is  ±  20%.  For  the 
accuracy  of  the  expansion  energy  at  the  diffuser  inlet  follows,  then 
a  figure  of  about  +  1.5%  or  for  the  expansion  Mach  number,  which  is 
about  a  square  root  function  of  the  flow  energy,  an  accuracy  of  somewhat 
less  than  ±  1%  follows.  Roughly  the  same  accuracy  applies  to  the  normal 
shock  recovery  ratio,  the  standard  of  comparison  for  the  diffuser  tests. 
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SECTION  VI 

RADIAL  DIFFUSER  DESIGN  CONSIDERATIONS 
1.  THE  MULTI-CHANNEL  RADIAL  DIFFUSER 

The  flow  problem  encountered  here  of  diffusing  a  supersonic,  radial 
flow  coming  from  a  cyl indrical ly  shaped  source  with  a  length  to  diameter 
ratio  near  one  or  larger,  is  without  a  commonly  known  precedence. 

Known  design  considerations  for  supersonic  diffusers  have  only  a  limited 
meaning  for  the  radial  diffuser.  For  example,  providing  a  constant  flow 
cross  section  for  the  radial  flow  by  proper  contouring  of  the  side  walls 
of  the  radial  flow  system  not  only  leads  to  complex  channel  geometries 
but  also  to  intolerable  flow  conditions,  since  at  desirable  L/D  ratios 
of  the  flow  source,  strong  axial  flow  deflections  become  necessary  near 
the  walls.  Subdivision  of  the  flow  becomes  an  absolute  necessity  to 
arrive  at  reasonable  flow  conditions  for  the  diffusion  process.  The 
subdivision  may  either  be  in  axial,  peripheral,  or  both  directions. 
Subdivision  i  i  only  one  direction  leads  to  undesirable  high  aspect  ratios 
of  the  diffuser  channel  cross  section.  The  present  design  provides 
subdivision  in  both  directions  resulting  in  a  multi-channel  diffuser 
(egg-crate  diffuser),  where  each  channel  can  be  given  a  cross  section 
shape  with  an  aspect  ratio  not  larger  than  two.  The  subdivision  into 
many  small  channels  allows  also  a  considerable  reduction  of  the  radial 
diffuser  diameter  since  the  supersonic  diffusion  process  depends  only 
on  a  given  length  to  diameter  ratio  of  the  individual  channel. 

The  physical  justification  for  applying  an  extreme  of  subdivision 
is  the  consideration  that  wall  friction  losses  play  only  a  minor  role  in 
the  supersonic  diffusion  process,  which  operates  on  the  principle  of  the 
so-called  pseudo  shock  consisting  of  a  system  of  supersonic  shocks 
recompressi ng  the  flow  very  nearly  in  the  same  degree  as  a  normal 
supersonic  shock.  In  this  shock  process  which  extends  for  a  length  of 
about  ten  diffuser  channel  diameters,  the  channel  walls  are  almost 
exclusively  exposed  to  subsonic  flow  conditions  only,  reducing  wall 
friction  to  an  almost  negligible  amount.  Thus  the  introduction  of  wall 
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surfaces  required  for  subdividing  the  flow  is  no  deterrent  to  an  effective 

diffusion  process,  as  long  as  the  Reynolds  number  of  the  channel  is 

4 

maintained  high  enough  (above  10  ).  The  subdivision  has  functional 
limitations.  For  instance,  if  the  diffuser  must  be  cooled,  as  it  is 
generally  necessary  for  a  high  energy  laser,  it  may  not  be  possible  to 
effectively  cool  the  wall  surfaces  in  case  of  a  high  degree  of 
subdivision.  Such  considerations  are,  however,  ignored  in  the  present 
investigations  in  preference  to  the  general  exploration  of  the  behavior  of 
a  multi-channel  radial  diffuser. 

2.  SPECIFIC  DIFFUSER  DESIGN  CONSIDERATIONS 

The  present  design  not  only  features  subdivision  in  two  directions 
but  also  in  two  steps.  The  flow  is  first  divided  axially  by  parallel 
plates  with  beveled  inlet  edges  as  seen  in  the  side  views  in  Figure  10 
(Section  A-A).  At  the  downstream  end  of  these  inlet  edges,  the  flow  is 
divided  also  in  peripheral  direction  by  the  pie  shaped  blades  seen  in 
the  front  view  of  Figure  10.  This  two-step  subdivision  is  dictated 
by  aerodynamic  requi rements .  The  first  division  provides  the  channel 
contraction  normally  desired  for  supersonic  diffusion.  Since,  in  this 
contraction  region,  the  flow  remains  supersonic,  wall  friction  has  an 
essential  influence  and  wall  surfaces  must  be  kept  to  a  minimum.  The 
shape  of  the  channel  cross  section  is,  on  the  other  hand,  not  important 
and  can  have  a  high  aspect  ratio  as  is  the  case  at  the  inlet  to  the 
diffuser  plates.  Where  the  contraction  ends  and  the  pseudo  shock 
systems  begins,  the  wall  friction  becomes  unimportant,  as  pointed  out 
before;  but  the  cross  section  shape  becomes  influential  and  should  be 
near  axisymmetric.  These  conditions  are  fulfilled  by  the  second  division. 
This  design  procedure  had  been  previously  applied  to  straight  flow 
supersonic  diffusers  in  Reference  3. 

The  Multi-channel  arrangement  allows  one  to  reduce  the  treatment  of 
the  radial  flow  diffusion  to  the  common  problem  of  the  straight  diffuser 
with  the  exception  of  the  inlet  flow.  The  inlet  edges  of  the  diffuser 
plates  are  circular  arcs  and  the  supersonic  flow  around  these  edges  must 
have  a  pronounced  three-dimensional  character.  The  flow  area  contraction 
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in  the  inlet  region  also  has  an  uncommon  shape  due  to  the  radial 
divergence  of  the  flow.  Figure  11  shows  the  area  reduction  at  the 
diffuser  plate  inlet  plotted  over  the  plate  radius  related  to  the  inlet 
radius  for  the  simple  wedge-shaped  inlet  edges  used  in  the  present  design. 
The  resulting  contour  nearly  has  the  shape  of  an  inverted  supersonic 
nozzle  for  parallel  exit  flow.  The  utility  of  this  shape  for  the 
diffusion  process  must  be  found  from  the  experiments.  No  attempt  has 
been  made  to  analyze  these  unique  inlet  flow  conditions.  Figure  11  also 
shows  the  diffuser  inlet  contour  for  the  theoretical  limit  of  inlet 
contraction  (M  =  3.514)  for  the  condition  that  the  cross-sectional  area 
decreases  linearly  with  the  radius. 

3.  THE  STRAIGHT  FLOW  SUPERSONIC  DIFFUSER 

In  the  following,  the  general  design  of  the  straight  channel  shock 
diffuser  is  discussed.  The  flow  area  schedule  of  the  diffuser  channel, 
i.e.,  the  diffuser  section  downstream  of  the  throat  of  a  supersonic 
diffuser  is  determined  by  two  factors:  flow  stability  and  performance. 

The  stability  requires  that  it  is  possible  to  position  the  beginning  of 
the  shock  system  very  close  to  the  inlet  of  the  channel  without  the 
shock  system  being  pushed  out  of  the  diffuser  by  a  small  disturbance. 

The  channel  flow  is  inherently  stable  if  an  upstream  movement  of  the 
shock  system  increases  the  pressure  recovery  of  the  diffuser.  In  this 
respect  a  constant  area  duct  is  stable  with  the  help  of  wall  friction. 

If  by  increasing  the  back  pressure  the  shock  system  moves  toward  the 
inlet,  less  supersonic  flow  is  exposed  to  wall  friction  and  the  pressure 
recovery  improves.  Thus,  the  shock  system  can  be  moved  upstream  until 
the  improved  recovery  balances  the  new  back  pressure.  Only  a  new  increase 
in  back  pressure  would  further  move  the  shock  system  upstream  and  if 
the  increase  in  back  pressure  persists,  the  shock  system  would  be  finally 
pushed  out  of  the  diffuser  channel. 

The  flow  stability  of  a  diffuser  can  be  increased  by  giving  it  a 
diverging  cross  section.  However,  divergence  reduces  also  the  overall 
diffuser  performance.  Thus,  some  compromise  is  necessary  to  find 
optimum  diffuser  conditions.  In  the  following,  a  graphic  method  based 
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on  equations  previously  derived  will  be  shown  and  this  method  allows  a 
check  on  the  stability  and  performance  of  a  diffuser  relative  to  its 
geometry  in  terms  of  its  inlet  to  exit  cross-sectional  area  ratio. 


In  Section  III,  the  equations  for  relating  the  flow  conditions  in  two 
cross  sections  of  a  non-cyl indrical  duct  with  wall  friction  have  been 
derived.  These  equations  can  be  conveniently  used  to  graphically  present 
the  exit  condition  of  a  shock  diffuser  in  relation  to  the  inlet  condition, 
the  duct  geometry,  and  the  wall  losses.  Equation  45  can  be  written  with 
the  help  of  Equation  35: 


1  1 

2  2 

M,0  +  •?  -  M2)  M-j  ( 1  +  ) 


'1 


?  ^  f ^ 

y  m2  +  1 T 


Y  +  i  •  t 


(79) 


In  this  equation,  the  functional  form  of  the  left  side,  which  is  given  by 
the  exit  conditions,  is  the  same  as  that  of  the  right  side,  which  is  a 
function  of  the  inlet  conditions.  Thus  a  graphic  presentation  of  the 
general  function 

M ( 1  +  M2)2 
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for  the  exit  gives  a  set  of  two  curves  which  correlate  inlet  and  exit 
Mach  number  for  a  duct  because  the  value  of  E  must  be  the  same  for  the 
inlet  and  the  exit.  Figures  12  through  14  give  examples  for  applying 
this  process  to  a  supersonic  diffuser  with  a  diverging,  a  constant 
area,  and  a  converging  channel.  Assuming  the  inlet  Mach  number  as  given, 
the  E-value  of  the  diffusion  process  is  found  from  the  "In"-curve. 

On  the  "Ex"-curve,  the  same  E-value  determines  the  exit  Mach  number  at 
the  subsonic  branch,  since  the  flow  underwent  complete  shock  diffusion. 

If  no  shock  diffusion  would  have  taken  place,  the  exit  Mach  number  would 
have  been  given  by  the  supersonic  branch  of  the  Ex-curve.  The  difference 
between  inlet  and  exit  Mach  number  in  the  absence  of  shock  diffusion  is 
then  caused  by  wall  friction  and  flow  divergence. 


The  ratio  between  inlet  and  exit  pressure  can  also  be  determined 
graphically.  Equation  36,  which  gives  this  pressure  ratio,  can  be 
written  with  the  help  of  Equation  35: 


(81) 


or  in  general  form: 
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The  coefficients  are  identical  with  those  used  for  plotting  the 
E-curves.  Thus  any  point  on  an  E-curve  can  be  assigned  a  value  given  by 

b  =  -Y-.:  •  c  +  l  (83) 

pa 

Bp-values  have  been  marked  on  the  E-curves  in  Figures  12  to  14. 

In  accordance  with  Equation  82  the  ratio  between  two  B  -values 

P 

belonging  to  two  different  E-curves,  but  having  the  same  E-value, 
provides  the  pressure  ratio  for  the  flow  between  the  two  cross  sections 
represented  by  the  respective  E-curves. 

To  find,  for  instance,  the  pressure  ratio  for  our  example  in  Figure  12. 

we  read  a  value  B  =  16.7  at  the  intersection  of  the  constant  E  line  with 
P 

the  In-curve.  On  the  subsonic  branch  of  the  In-curve,  we  find  =1.26.  The 

ratio  B  /B  =  13.25  is  then  the  static  pressure  ratio  between  inlet  and 
P  P 

exit  of  our  diverging  diffuser  duct. 

To  check  on  the  stability  of  the  flow,  i.e.,  whether  the  shock  system 
moves  uncontrollably  upstream  by  an  increase  in  back  pressure,  we 
investigate  the  two  cases  where  the  shock  system  starts  at  the  channel 
inlet  and  where  it  starts  further  downstream.  In  the  second  case,  we 
have  supersonic  flow  extending  over  a  portion  of  the  diffuser  channel. 

The  pressure  ratio  across  this  channel  portion  is  obtained  from  the  Bp  values 
at  the  intersections  of  the  constant  E  line  with  the  supersonic  branches 
of  the  two  E-curves 

B  =  16.7  and  B  =  16.0  then  Pr  / PT  =  16.7/16.0  =  1.044 
p  p  Ex  In 

The  inlet  Mach  number  to  the  downstream  portion  of  the  diverging 
diffuser  channel  where  the  shock  system  is  now  located  is  the  exit  Mach 
number  of  the  upstream  portion  where  the  flow  is  supersonic.  In  Figure  12, 
we  also  obtain  with  this  new  inlet  Mach  number  a  new  E-value  on  the 
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In-curve  and  the  exit  Mach  number  after  complete  shock  diffusion  from 
the  Ex-curve.  From  the  pertinent  B^-values  (13.1  and  1.28)  we  obtain 
the  pressure  ratio  across  the  shock  system 


Pr  (B  )r 
Ex  _  '  p'Ex  _ 

PT  ‘  (B  )' 

In  p  In 


10.23 


The  overall  pressure  ratio  across  the  diffuser  channel  is  then 
10,23  *  1.044  =  10.68,  We  see  that  this  pressure  ratio  is  lower  than 
that  resulting  for  the  shock  system  located  upstream.  Since  the  pressure 
recovery  ratio  increases  with  the  shock  system  moving  upstream  the 
diffuser  is  stable. 


Table  3  lists  all  magnitudes  of  interest  derived  from  Figures  12 
through  14  for  the  diverging,  constant  area,  and  converging  case.  For 
the  converging  channel  in  particular,  we  find  that  the  pressure  ratio 
for  the  shock  positioned  downstream  is  now  higher  than  when  the  shock 
system  begins  at  the  channel  inlet;  i.e.,  the  flow  in  this  converging 
channel  is  unstable. 

The  above  considerations  about  diffuser  area  scheduling  provide  the 
following  guidelines  for  shaping  the  diffuser. 

1.  The  diffuser  should,  in  the  common  way,  begin  with  a  converging 
section  in  which  kinetic  energy  is  converted  to  pressure  energy  with  a 
minimum  of  shock  losses.  The  throat  of  the  diffuser  is  at  the  end  of 
the  converging  section. 


2.  A  diverging  section  should  follow  the  converging  section  for 
stabilizing  the  shock  front  as  close  to  the  throat  as  possible.  Its 
length  should  not  be  greater  than  necessary  for  stabilizing  the  flow 
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3.  For  the  remaining  length  of  the  shock  system,  the  diffuser  duct 
should  be  converging  again  or  be  at  least  of  constant  area. 

NOTE:  A  supersonic  diffuser  commonly  consists  of  two  sections, 
a  converging  inlet  section  and  a  diverging  or  constant  area  section 
for  the  shock  diffusion  process. 

Experimental  experience  can  only  show  the  benefits  of  dividing  a 
supersonic  diffuser  into  three  sections  as  outlined  above.  Figure  10 
gave  the  design  of  a  radial  diffuser  plate  for  the  present  diffuser  test 
In  this  figure,  two  versions  are  indicated;  one  where  the  throat  is 
followed  by  a  constant  area  section  and  one  where  one  wall  of  the  fol- 
lowing  section  diverges  by  about  1°.  Additional  adjustments  of  the 
diffuser  channel  geometry  are  to  be  accomplished  by  epoxy  layer  built  up 
on  the  channel  walls  of  the  machined  diffuser  plates. 
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SECTION  VII 
CONCLUSION 

Supersonic  radial  outflow  expansion  systems  are  afflicted  with 
geometric  constraints  which  make  their  design  and  also  their  operation 
more  complex  than  is  common  for  straight  flow  systems.  To  solve  this 
design  problem,  design  charts  have  been  developed  which  relate  the 
principal  geometry  of  ring  nozzles,  the  basic  elements  of  the  expansion 
system, to  the  flow  conditions  in  terms  of  the  isentropic  expansion  Mach 
number,  weight  flow  rate,  upstream  and  downstream  pressures,  and  throat 
Reynolds  number.  Design  charts  of  this  kind  have  been  prepared  for  four 
different  gases  (He,  air,  CO^ ,  and  Freon  12).  The  chart  for  air  has  been 
used  to  design  two  radial  flow  nozzle  systems  (M  =  2.3  and  3.3)  for  use 
in  a  test  rig  to  explore  radial  flow  diffusers. 

Radial  flow  diffusers  used  for  laser  systems  feature  unusually  large 
L/D  ratios.  This  geometric  condition  affects  the  operation  of  the 
diffuser  cest  rig  considered  here  twofold.  It  requires  subdivision  of 
the  flow  into  many  small  expansion  units  and  prevents  the  access  to  the 
expansion  system  exit  for  surveying  the  flow  by  means  of  probes. 

Subdivision  introduces  flow  losses  and  simple  potential  flow  cons i dera tions 
are  no  longer  adequate  for  determining  the  expansion  performance.  Since 
probing  of  the  flow  is  also  not  possible,  the  expansion  exit  conditions, 
which  are  the  inlet  condition  of  the  diffuser,  must  be  determined 
analytically  with  careful  consideration  of  all  possible  flow  losses. 

The  analytic  process  which  accounts  for  nozzle  wall  friction,  losses 
due  to  sudden  flow  expansion  on  the  nozzle  edges,  and  the  radial 
divergence  of  the  flow  is  developed  here  for  a  ring  nozzle  system. 

For  the  specific  nozzle  design  to  be  used  for  the  diffuser  test  rig,  the 
analysis  finds  an  effective  Mach  number  of  M  =  3.514  (for  isentropic 
expansion  M  =  3.688).  For  evaluating  diffuser  test  results,  the  measured 
pressure  recovery  is  compared  to  the  normal  shock  recovery  pertaining 
to  the  effective  Mach  number  indicated  above. 
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For  an  effective  diffusion  of  the  flow  in  a  radial  flow  system, 
subdivision  of  the  flow  field  into  small  units  is  necessary.  Providing 
a  near  constant  area  flow  path  as  required  for  supersonic  shock  diffusion 
leads  to  complex  flow  fields  in  case  of  a  radial  flow  field  with  no 
subdivision.  This  is  particularly  true  for  the  large  L/D  ratio  systems 
needed  for  radial  flow  lasers.  Since  the  walls  are  exposed  only  to 
subsonic  velocities  in  an  effective  shock  diffuser,  subdivision  poses 
no  flow  loss  problem  as  long  as  the  channel  Reynolds  number  does  not 
become  too  low.  Subdivision  must  be  minimized  ahead  of  the  diffuser 
throat,  where  velocities  are  still  supersonic  on  the  wall  (subdivision 
should  be  also  minimized  for  any  suosonic  diffusion  following  the  shock 
diffusion  since  there  is  no  longer  any  inherent  need  for  subdivision 
in  a  radial  subsonic  diffuser).  Subdivision  can  become  an  operational 
problem  for  laser  applications  where  it  may  interfere  with  cooling  the 
di f fuser . 

Subdivision  of  the  radial  diffuser  into  many  small  units  makes  it 
possible  to  apply  a  standard  geometry  to  the  individual  channels.  Except 
for  the  diffuser  inlet  where  a  complex  three-dimensional  flow  field 
exists,  common  design  considerations  apply  to  the  subdivided  radial 
diffuser.  Flow  stability,  which  is  the  principal  criterion  for  the 
shock  diffuser  design,  is  established  in  a  duct  if  the  pressure  reco\ery 
across  the  shock  system  increases  when  the  shock  system  is  moved 
upstream  toward  the  diffuser  throat  by  increasing  the  diffuser  back 
pressure.  In  this  respect  wall  friction  makes  the  flow  stable  in  a 
constant  area  diffuser  channel.  Stability  is  increased  by  providing  a 
diverging  channel  at  the  expense  of  the  diffuser  performance.  A  possible 
compromise  is  to  provide  a  short  diversion  after  the  diffuser  throat  and 
continue  with  a  constant  area  channel. 
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Figure  1 .  Schemp 
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Figure  2.  Ring  Nozzle  Design  Chart  for  Helium.  Dimensions: 
P(Newton/m2) ;  H(m);  r(m);  T(CK);  W(kg/sec) 


/ 


AFWAL-TR-80- 3028 


P  ■  H 


o 


Figure  3.  Ring  Nozzle  Design  Chart  for  Air.  Dimensions 
P(Newton/m^) ;  H(m);  r(m);  T(°K);  W(kg/sec) 


Figure  4.  Ring  Nozzle  Design  Chart  for  CC^-  Dimensions 
P(Newton/m^) ;  H(m);  r(m);  T(°K);  W(kg/sec) 
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Figure  5.  Ring  Nozzle  Design  Chart  for  Freon  12.  Dimensions: 
P(Newton/m^) ;  H(m);  r(m);  T(°K);  W(kg/sec) 
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Figure  6.  Design  of  a  Ring  Nozzle  Expansion  System  for  M 
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Figure  7.  Scheme  of  a  Radial  Outflow  System 
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Figure  9.  Distribution  of  the  Wall  Friction  Forces  Along  a 
Supersonic  Expansion  Nozzle  for  a  Constant  Hall 
Friction  Coefficient 
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Figure  11.  Wall  Contour  and  Cross-Sectional  Area  for  Two  Different 
Radial  Diffuser  Inlet  Geometries.  (Wedge  shaped  inlet 

{ - )  is  for  a  throat  contraction  ratio  0.86;  constant 

rate  of  inlet  area  decrease  ( - )  is  for  a  contraction 

ratio  0.69 


